ASYMPTOTIC CORRELATIONS FOR GAUSSIAN AND WISHART 
MATRICES WITH EXTERNAL SOURCE 



o 
o 



> 

(N 

o 
o 
o 

-)— > 



X 



PATRICK DESROSIERS AND PETER J. FORRESTER 



Abstract. We consider ensembles of Gaussian (Hermite) and Wishart (Laguerre) N x N her- 
mitian matrices. We study the effect of finite rank perturbations of these ensembles by a source 
term. The rank r of the perturbation corresponds to the number of non-null eigenvalues of the 
source matrix. In the perturbed ensembles, the correlation functions can be written in terms 
of kernels. We show that for all A'', the difference between the perturbed and the unperturbed 
kernels is a degenerate kernel of size r which depends on multiple Hermite or Laguerre functions. 
We also compute asymptotic formulas for the multiple Laguerre functions kernels in terms mul- 
tiple Bessel (resp. Airy) functions. This leads to the large N limiting kernels at the hard (resp. 
soft) edge of the spectrum of the perturbed Laguerre ensemble. Similar results are obtained in 
the Hermite case. 
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2 patrick desrosiers and peter j. forrester 

1. Introduction 

We study two specific ensembles of N x N random hermitian matrices X in which the prob- 
abihty distribution function (p.d.f.) takes the general form 

P(X|A) OC e-T'^^W-TrAX^ (1_1) 

where A is a given N x N hermitian matrices and V^(X) is a matrix-valued function of X. These 
matrix ensembles have been considered by many authors, notably Brezin and Hikami [S], Guhr 
and Wettig |2n|l21j. Zinn- Justin '33 ; and more recently, Baik, Ben Arous and Peche [Hl3ll29]. 
El Karoui ^2j, Bleher and Kuijlaars[^^, Imamura and Sasamoto |23] . 

Using the celebrated Harish-Chandra |22] (or Itzykson-Zuber [23j) formula, 

/ e^("^"'«'(du) = nr(i) /■^"^•'^-ly;;-:; (1.2) 

J ._i det a^ ij=i,...,Afdet 6^ \i,j=i,...,N 

where (dXJ) stands for the normalized Haar measure, one readily shows that the eigenvalue joint 
p.d.f. is 

P^(x|a) = ^det[e-^-^]f,_^n^"''^'"^ U ^^^' (l-^) 

where x = [xi, . . . , xn) denotes the random eigenvalues of X while a = (ai, . . . , ajy) denotes the 
fixed eigenvalues of A. The quantity Zjv is the normalization. We suppose that xi > . . . > xn 
and that Xj £ J CI R. An essential feature of ()1.3() is that it is of the form 

P^(x|a) = ^ n —^ det Uxk)]f,^, det [Vj{xk)]fk=i (1-4) 

^ l<J<k<N "^' "j' 

(set ^j{x) = x^~^e~°'^^~^^^' , r]j{x) = x^~^, and make use of the Vandermonde determinant 
evaluation in the latter). We are thus dealing with determinantal point processes, which can be 
studied as biorthogonal ensembles in the sense of Borodin |H] . 

In this article, we consider the large A^ behavior of the correlations in two biorthogonal 
ensembles having an eigenvalue p.d.f. of the form H1.3() . 

1) V{x) = x^ J = (— cxD,oo) Gaussian (Hermite) Unitary Ensemble; , . 

2) V{x) = X — aln(x) J = (0, oo) Wishart (Laguerre) Unitary Ensemble. ' 

To simplify subsequent integral representations, a will be taken to be a non-negative integer; 
however we expect our final formulas, appropriately interpreted, to be valid for general real 
a > -1. 

Both the choices in (|1.5|) can be realized by matrices with Gaussian entries. Consider case 1). 
Let A be an A^ X A^ Hermitian matrix with distinct eigenvalues ai, . . . , a^- Let Y be an element 
of the GUE, and thus be an A^ x A^ complex Hermitian matrix with probability density function 
proportional to e~ . With t a positive scalar, it is a well known consequence of p.2|) (see 
e.g. jl2j) that the eigenvalue p.d.f. of 

-A + VtY (1.6) 

is equal to 

izl);!5^!:!^det[e-(^'-^)^/-i--, ,^i^^[M:i-wv .,. 

m2nt)m ''*^' ^^•^^"■■■''^detK-^].,.!,...,^- ^ ^ ^ 

Simple manipulation and use of the Vandermonde determinant identity reduces (|1.7|) in the case 
t = 1 to dUHl) with V{x) = x^. 
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Now consider case 2). Let W be a n x A^ (n > A^) complex Gaussian matrix with real and 
imaginary parts having variance 1/2, mean 0, and let B be an A^ x A^ positive definite Hermitian 
matrix with distinct eigenvalues bi, . . . ,bN. Again it is a straightforward consequence of (|1.2|) 
d Cni EOl that the eigenvalue p.d.f. of 

B-V2wtwB"i/2 (1.8) 

is equal to 

(-l)«(A'-i)/2 IL 6f A?-" „ (A. - A,) ,., 

W<. nin-iv + i-D-JI^^W^j^-J^'l^ '!«- " ('■») 

and thus with 6j = 1 + aj, Xj = Xj, is proportional to ()1.3|) with V{x) = x — {n — N) logx. 

A specialization of (|1.8() of interest in mathematical statistics [2S] is to take B = diag(6i, . . . ,br, 
1, . . . , 1). Then the matrix X := Wl^B"^'^ has columns j, (j = 1, . . . , r) with variances 1/2-^/6^, 
while the entries in the remaining columns all have variance 1/2. Regarding X as a data matrix, 
this corresponds to so called spiked data. It was shown in [^ that with the scaled variables Sj 
and Xj specified by 

''^ = ^ + 2{2N)y^ '^ (i = l>--->0 x,=4A^ + 2(2Ar)i/3x, (j = 1, . . . , fc) (1.10) 

the scaled fc-point correlation 

1 s h 

(1.11) 



li,u=l, 



( ^ \^ 


(xi,...,Xfe) 


V2(2Ar)i/3J ''w 


tends in the hmit N ^ oo, n/N ^ 1 to 


p;;*^^^>(Xi,...,X,):=det 


1C'°'\X^,X,) 



where 



with 



(1.12) 

J ^,u=l,...,K: 



IC'°^^{X,Y) =K'°^\X,Y) + ^Ai{X)Ai^'\Y) (1.13) 



i=l 



,,-„«(^,,,),Ai(A-)Ar(y-)-Ai'(X)Ai(r)_ i^^^j 

A^ — y 

Ai^^x) = /■ l^e— +^V3pj(, _ ,^)^i, (1.15) 

•^{ai,...,Sj} 

Ai«(a;) = (-1)^ / ^^e-^^+^'/^lliv + .,). (1.16) 

See Fig. [21for the definition of the contours ^|si,...,sj and A. (The same limit ()1.12|) was found 
for N ^ oo with n/N ^ 7^, 7^ G M"*", provided the scales in (|1.1()|) are suitably modified.) The 
significance of the scale p.l()|) is that it centres the coordinates about the neighbourhood of the 
largest eigenvalue, and scales these eigenvalues so that their spacing is of order unity. Because 
the eigenvalue density is not strictly zero beyond the expected position of the largest eigenvalue, 
this is referred to as the soft edge, and thus the reason for the superscript 'soft' in (|1.12j) and 
(dlHl). 
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Subsequent to the derivation of (|1.12)) as a limiting correlation at the soft edge for matrices 
(|1.8j) with B a rank r perturbation of the identity, it was shown in [22] that (|1.12)) is also 
the limiting correlation at the soft edge for matrices (|1.6jl with A a rank r perturbation of 
the zero matrix. Explicitly, if in (|1.6|) we take i = 1 and A = diag(ai, . . . ,aj.,0, . . . ,0), with 
ai = • • • = Or = 1, and introduce the scaled variables Xj = \/2N^'^{Xj — \/2N) {j = 1, . . . , k), 
then in the limit N ^ oo the scaled /c-point correlation 

(x/2ivV6)V(fc)(xi,...,Xfc) 

tends to H1.12() with si = • • • = s^ = 0. 

In both the studies |5, ^\ the structure H1.13() appears after the asymptotic analysis of the 
respective double contour representation of the general finite N, (fi, v) element in p.lH|) (which 
is called the correlation kernel). It is an objective of this to identity the structure analogous to 
(|1.13|) in the finite A^ p.d.f. s (|1.7|) . (|1.9|) . Moreover for (|1.9|) we identify a hard edge scaling (a 
scaling of the smallest eigenvalues in (|1.9j) ) to the limiting A;-point correlation 



where 



hard{h,}^^^^ . . . ,X,.) := det f/C^-'i(X^,X,) (1.17) 



/i,!/ = l, 



(a+r)/2 



with 






^a {^^^) — 2(x - Y) ' u-J^yj 

J«(x)= / —— , (1.20) 

„.,,,, ^ / ju,e--/^--na(^--%), ,1,21) 

^ ^ 7 27ri «;"+'• ^ ' 

C{o} 

Here Csq^^^ ^^a is a simple closed anticlockwise contour encircling the points hi and the origin, 
while Cjo} is a simple closed anticlockwise contour encircling the origin (see Fig. ^ . 

The correlation functions obtained here are multiparameter generalizations of the correlations 
in the Hermite and Laguerre unitary ensembles. The latter quantities are reviewed in the next 
section. In §3, we present some elements of the theory of multiple orthogonal functions, which 
provides the natural setting for studying the matrix ensembles with p.d.f. of the form (|1.3|) . We 
show in §4 that /c-point correlation function can be written as det[i^(xi,Xj)]jj=i^...^fc for a given 
kernel K(xi,Xj). The finite rank, r say, perturbations of the Hermite and Laguerre ensembles 
are considered in §5. We show that even for matrix ensembles of finite size A^ the kernel takes the 
form exhibited by (|1.13|) . (|1.18|) . Explicitly K]\f{x,y\a) = K]^-rix,y\0) + J2'i=i /i(^)/i(j/)) where 
fi and fi denote multiple functions that generalize the Hermite or the Laguerre polynomials. In 
§6, we consider the large A^ limit and obtain the scaled kernels (|1.13|) . (|1.18j) . This is achieved 
by showing that the multiple Laguerre functions tend to multiple Bessel functions at the hard 
edge, while both the multiple Laguerre and Hermite functions tend to the same multiple Airy 
functions at the soft edge. The article ends with a short discussion. 
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2. Null case 

The null case corresponds to ensembles of random matrices with unitary symmetry (see for 
instance ^lE^j). It is defined by 

a^0 = {0,...,0} (2.1) 

in (|1.3|) . In such a situation, ni<i(%' ~^j) det [e~'^''^']f^ i tends to ni<?(^i ~^i) ' ^'^d the fc-point 
correlation function can be written as a fc x /c determinant. 

In order to be more concrete, let us introduce Pn{x) = x" + c„_ix"~^ + . . ., an orthogonal 
polynomial with respect to the weight w{x) := e~'^^^' supported on the real interval J; that is, 

dxw{x)pj{x)pk{x) = WpjW'^dj^k- (2.2) 

J 
Let also (Christoffel-Darboux formula j31j ) 

Af-l 



J^ / ^ / / X / N V^ Pn{x)Pn{y) 

KN{x,y) := yjw{x)w{y) 2^ 



II 119 

n WJr'nW 

n=0 



^ ^yw{x)w{y) PN{x)pN~l{y) -pN-l{x)pN{y) 

\\PN-i\? x-y 

be the reproducing kernel on the interval J. It satisfies 



(2.3) 



dyKN{x,y)KN{y,z) = Kn{x,z), I dxKN{x,x) = N (2.4) 

J J 

and 

1 ^ 
PN{y^) = ^W ^(^i) n (^^ ~ ^j)^ " ^^^ \KN{xi,Xj)]^^.^^ . (2.5) 

^ j=i i<i<fc<A' 

These properties together with the Laplace expansion of the determinant imply that the A;-point 

correlation function can be expressed as a k x k determinant of the kernel, 

Nl f f 

Pik){xi,---,Xk) ■■= , / dxk+i--- / dXNPN{xi,...,XN) 

J J 

= det [KN{xi,Xj)]^.^-^ (2.6) 

For the Hermite and Laguerre ensembles, the monic polynomials of degree n can respectively 
be expressed in terms of the Hermite and the Laguerre polynomials, 

1) Pnix)=2-^Hnix) |K||2 = ^2-"n! if J = (-00,00), , . 

2) p„(2;) = (-l)"n!L^(x) ||p„||2 = n!(n + a)! if J=(0,oo). ^^^ 

In the Hermite case, we thus have 

^-i.^W)/2 Hr,{x)HN-i{y) - HN-i{x)HN{y) ,_, 

^^^^' y^ = WvW^Ty. l^y ' ('•') 

while we have, in the Laguerre case, 
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Of particular interest are the scaled kernels ^31; relating to the edge of the spectrum, which 
are defined by 

^scaied^^^ y) := lim B{N)Kn {A{N) + B{N)X, A{N) + B{N)Y) , (2.10) 

N~*oo 

where 'scaled' stands for 'hard' or 'soft'. When the functions A and B are suitably chosen, one 
gets kernels that are independent of A^ and that can be expressed in terms of simple functions. 
In the Hermite ensemble, there is only one distinct edge scaling. 



X = A{N) + B{N)X = Vm + X/V2N^/^, softedge. (2.11) 

Two edge scaling are possible for the Laguerre ensemble, 

X = A{N) + B{N)X = i''/^''\ ,,,3 ^^^^^^S^' (2.12) 

^ ^ ^ ^ \4N + 2{2Ny/^X, softedge. ^ ' 

By substituting the corresponding scaling of Eqs. ()2.11|) and 1)2.12(1 into Eq. ()2.1U|) the scaled 
kernels (fm)) . (fTT^ result. 

3. Multiple biorthogonal polynomials 

3.1. Generalities. We introduce the multiple Hermite and Laguerre polynomials mSl[7j. Note 
that we use slightly different notations and definitions than those of the latter references in order 
to simplify the comparison between our asymptotic formulas and the well known formulas as 
revised above in the null case. 

Let us suppose that within the set 

a = {ai,...,aAr}, (3.1) 

only D eigenvalues are distinct. Then we may write 

a = b-:={6-\...,6™-}, (3.2) 

which means there are rrii eigenvalues in a that are equal to 6j (i = 1, . . . ,D). Let also |m| := 
X^i=i "^« — ^i that is, m is a composition of non-negative integers, of weight N and of fixed 
length D. 

The ith multiple orthogonal polynomials of type I associated to a given set of parameters 
b, denoted by Qm,iix) or Qm,i{x\h) (but usually, the dependence on b is kept implicit), is a 
polynomial in x of degree rrii — 1. It is used to build the function Qm{x) := Yli=i ^~^^^Qm,i{x) 
such that 

dxw{x)x^Q mix) = < ' ' ' ' (3-3) 

[1, J = |m| - 1. 

The multiple orthogonal polynomial of type II and of degree |m| is denoted by Pm{x) or Pin(x|b). 
It satisfies 



/ 



dxw{x)e-'''''x^P,rvix) = 0, j = 0,...,mi-l, (3.4) 



J 
where i = 1, . . . ,D. Combining the two previous equations, we see that 



dxw{x)Pni{x)Qnix) 

J 



0, |m| < |n| — 1, 

1, |m| = |n| — 1, (3-5) 
0, rrii > rii — 1 Vi. 
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In other words, for a given set of parameters b and a given composition m of weight N, the 
polynomials of type I and II allow us to build a family of biorthogonal functions, 



dx w{x)pi{x)qj{x) = 5i^. 



for all z, j = 0, 



,N -I 



by choosing 

Pi{x) = P^^ (x) , Qiix) =Qf^^^^{x), 

and if for instance, the compositions /Xj are given by 

/ fij^ \ (mi 



(3.6) 



(3.7) 






1712 



ma 
mi — 1 ?7i2 m,3 

m,i — 1 ?7i2 — 1 m,3 
mi — 1 771-2 — 1 'iT-a — 1 











rriD \ 
mo 
mo 
mo 

J 



(3.8) 



where it is understood that mi > ?n-j_|_i. 

3.2. Multiple Laguerre functions. The multiple Laguerre function of type I is [3 cf. Eq. (3.10)] 

^ " dz e-^'^(l + z)l'"l+"-i 

1=1 



•^mV-^J ~ / ^ ^ ■'-'m,iV'' 



C, 



{-1} 



2^i nf=i( 



a S 



(3.9) 



where C^ = CJ^^-^ b | denotes a closed contour which encircles positively bi, . . . ,bD but not 
the point —1 (See Fig. ^. Note that encircling the point —1 is forbidden only when |ni| + a — 1 
is negative. This means in particular that for a > —1, the contour C^ ^ can be replaced by Cb. 
Obviously, the ith polynomial 'C^^(x) has degree m,j — 1. The multiple Laguerre polynomial of 
type II has degree \m\; it can be defined by [21 cf. Eq. (3.5)] 



C^{x) 



|m| + a)\ 



m ! 



dw 



27ri ■u;l'"l+°+i 



C{o} 



D 

m 

i=l 



W 



I -hi 



.rrii 



a G 



(3.10) 



where C/q} is a counterclockwise contour around the origin. Our definitions, modifying those of 
[7j are motivated in part by the first of the following results, while the key orthogonality type 
results remain essentially unchanged. 

Proposition 1. Let L"{x) denote the Laguerre polynomial of degree n. Set a a non-negative 
integer and bi, . . . ,bD > —1- Then the multiple Laguerre functions satisfy 



lim. CZ,(x) = L?L, i(x) and lim £J 

b^O |m| IV / j^^Q 

Moreover ,7.,, Thm. 3.2] 

'o, 



dxe-''x''+^C^{x) 







LU^)- 



j = 0, ...,|m| 



-l)l'"l-i(|m|+a-l)! 



T-D 



ur=ii^ + hy 



m 



1, 
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Figure 1. Contours C/„; and C 



{v} 



'W 



{-„j ana Cr'^^' ti > ^^ ^^^ complex z-plane. 

5(z) 













C^'^ • 


^) 














■,6d} 




r>^. 


V 


/: 




t 


u ^ 


^ 


62 ••• 


foy 


1 











3fi(z) 



and T, Thm. 3.1] 

00 

/"dxe-^-''»^x"+-'£^(x) = /or i = l,...,D and j = 0, . . . ,mi - 1. 



Proof. One easily proves the two first properties by considering the integral representation of 
the Laguerre polynomials [SH, 



L^(x) 



dw e-^'" ,, 
TT l + ^t' 



'{0} 



and the readily verified formula 

n!(-x)-"L;;-(x) = (n-a)!L-_Jx), 
both relations being valid for a G Z. Then, using Eq. (|3.9() . we have 

dz (l + z)l™l+"-i 



(3.11) 



(3.12) 



(ixe"^x"+^£^(x) 



Cb u 



dxe ^ ^^x"'^^ 



(a + j)! 



,{-1} 



dz (l + z)IH-i-' 



Note that, on the first line, we assumed bi > —1 and 3f?(z) > —1; this allowed us to do the 
2;-integration. For j = |m| — 1, we use the residue theorem and get 



dxe~^x"+J£^(x) = -(|m| + a - l)!Res^=_i 



D 



n(^ - ^^y 



.i=l 



as desired. For j < |m| — 1, we consider (1 + zy"^' ^ ^ ni=i(-2 ~ h) 
Thus, by virtue of the residue theorem. 



(|m| +a- 1)! 
► z~^~'^ when z — > 00. 



dxe ''x"+-'£^(x) oc Res2=oo [z ^ ^] = when - 1 < j < |m| - 1. 
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The same tricks are used to prove the second orthogonahty property. Exphcitly, 



I 



dxe-^-''»^a;"+^£"(x)oc 



dw Uf=i{w - I - bj 



27ri 



w 



|m|+Q:+l 



dx e-^-^'^+^^-x^' 



'{0} 



oc 



yi(w)<l+bi 

dw{w-l- bi)-^-' Uf=i{w - 1 - bjY 



/dw 
2^' 



w 



|m|+a+l 



'{0} 



~j~a-21 



OC Res^=oo [w ^ " ^\ = 0. 
when —l<j<mi and a > — 1, as proposed. 

Comparing with the general formahsm of §3.1 we have 

w(x) = x°'e~^ 



n 



Qin{x) 



(-i)|'"|-'nf=i(i+&^ 



(m + a- 1)! 



-C^{x) 



P^ix) = C^{x) 



(3.13) 



with corresponding biorthogonal functions specified by l|3.7|) . 

The multiple Laguerre functions also satisfy certain differential equations ^. Indeed, if 



d 



then 



Consequently, 



and 



Vt,J{x) := -e'^d^ (e-*V(x)) = {t - d^)f{x), with d, - ^^ , 



= V_a^,x---'D. 



CZ,{x) 



(_l)|m| 



m 



Jm|+Q 



(3.14) 

(3.15) 
(3.16) 

(3.17) 



When rriD parameters in a are equal to zero, aiml- 



niu 



'^|m| = say, after some 



simple manipulations we can show that the multiple Laguerre functions of type I and II can be 
respectively interpreted as anti-derivatives and derivatives of Laguerre polynomials. Explicitly, 
for r = |m| — mD, we see from the integral representations (|3.9|) . (|3.22|) and the differentiation 
formula ()3.15|) that 



|ni|— r- 



1 — ^—ai,x ' ' ' ^—ar,x •'-mV'^J 



(3.18) 



while the integral representations 1)3. 9|1 , ()3.22|) , the identity (|3.15|) and the differentiation formula 
H3.15() show that 



^"^S.(^) = (-ir ^'"!' /^' Pa.,. • • • I^a.,. X^+^L^^lix) 



m ! 



|m| — rV 



(3.19) 
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3.3. Multiple Hermite functions. The multiple Hermite polynomials of type I, written 'Hui,ii-c), 
are obtained from the function cf. (2.9)] 



D 



-Hmix) = Y,e-'^^n^4x) = (-2)l-l-i(|m| - 1)! 



4 = 1 



2^nf=i(^-^.r 



(3.20) 



where Cb is a simple contour which encircles the points {5i, . . . j^d} anticlockwise. The (non- 
monic) multiple Hermite polynomial of type II is given by [3 cf. (2.2)] 






(3.21) 



One can verify that the degrees of "^111,1(2^) and 7im{x) are respectively rrii — 1 and |ni|. 

The multiple Hermite functions (polynomials) can be obtained as a limit of their Laguerre 
counterpart. 



Proposition 2. Let X be real numbers such that |v2a-'^| < a with a S N. Then, as a ^^ 00, 
(2a)(i-l'"l)/2£^(Q, + V2^X) 



(_l)|mhl ^ / 1 

b.^bjV2^ = 2H-i(|m|-l)!^'"^^^ ^ ° V7^ 



and 



(2a)-l'"l/2£^(a + V2^X) 



f-Dlml / I 

= I I „ ^m(^) + O ^ 



Proof. The proof of the first asymptotic formula is the simplest. From the definition of the type 
I Laguerre functions we have 



(2a)(i-l'"l)/2£^(a + V2^X) 

bi^bi/V^ 

By changing z 1— > z/\/2x^, we get 

(2a)(i-l'"l)/2£^(a + V2^X) 



(2a) 



(l-|m|)/2 



dz e-°^-v^^^(l + z)l™l+"-i 



■'b/-/2a 



2^i nf=i(^-b^/V2^) 



6ii-^bi/\/2« 7 27ri 
Cb 



dz e-^^-V^^(l + z/^/2^)°+l'"l-i 



r-D 



nr=i(^-&o^ 



From the binomial development (1 + z/\/2a)"^l"*l ^ and the Maclaurin expansion of e V"/ ^^ 
we easily find that 

e-^/^^(l + z/^/2^)°+l'^l-i = e-^^^(l + z/^/2^)" + O (^^ 

when a ^ 00. This immediately implies the first asymptotic formula. 

The second asymptotic relation is established by steepest descents. The integral formula 
(|3.in|) shows 

£^(a + V2^X) = M + ^(« + V2^X)- / p,e-fMg{w), 
|ni|! J ZTTi 

C{o} 
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where 

D 

f{w) =w-lnw and g{w) = e-°Mi+\/57^^)eV^^«' JJ(ty - 1 - bi)"^' w-^"'^-\ 

We point out that aw dominates ^TlaXw by hypothesis. Thus, the function /(if) determines the 
principal contribution to the integral when a becomes large. fiyS) possesses one simple saddle 
point at Wo = 1 ; that is, f'{wo) = but f"{wo) = 1 7^ 0. The steepest descent contour is fixed 
by the condition (w — wq)'^ f" (wq) = minimum, which means arg(u;) = ib7r/2. We consequently 
choose the contour Cq of w such that the latter variable starts at 1 — ic (c is real and positive), 
passes through the point w = 1, goes to 1 + ic and finally returns to the starting point after 
encircling that origin. We now set v = \/^a{w — 1). When q — > 00 

af{w) =a + v'^/4 + ' 



a 

D 



{2a)\^\"g{w) = e^^+^' \{{v - V2^b,) + O (-^) ; 

and the steepest descent path in the complex w-plane becomes the interval (— icx),ioo) on the 
imaginary axis. We also know from Stirling's approximation that 

(|m| + a)! = ^/2^e-"al'"l+"+^/2 + q f ^ 



Therefore, we have proved 

ioo ^ 

bi^bi/V2a il""! m ! J ZTTl ^^ 

-ioo *=1 

which is the desired expression. D 

Proposition 3. Let Hn{x) he the Hermite polynomial of degree n. Then the multiple Hermite 
functions satisfy 

lim 7^111(2^) = H\^\_i(x) and lim Hmix) = H\^\(x). 



0, j = 0,...,|m|-2 

A/7r(|m| - 1)!, j = |m| - 1, 



Moreover [3 Thm 2.3] 

00 

/ dxe~^ x^Hm{x) 

—00 

and 7, Thm 2.1] 
00 

/ dx e"'^^"^*^x-'Wm(x) = for i = l,...,D and j = 0, . . . , mj - 1. 

—00 

Proof. The two first properties are straightforward consequences of well known integral repre- 
sentation of the Hermite polynomials j31j 



100 

2/ 



H^ix) = 2-nl I ^'^^^ = V^e^' [ ^^"e— +-^/l (3.22) 

^ ^ ' 27ri ■w'^+i ^ / 27ri ^ ^ 



c. 



{0} 



(l + W-xr = e-" +0 

*' a V A/a 
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We can prove the orthogonality relations by using Proposition El and substituting 

in the last equations of Proposition ^ D 

Suppose in (|3.2U|) and (13.21(1 that 6j = Oj, m,j = 1 (z = 1, . . . , r) while 6^+1 = 0, ra^j^x = |m| — r. 
Then we see from plKj) and (|3?22|) that 

dml - 1)! 

P_„,,, • • • V_ar,xnm{x) = {-^Y'jr^ ^^|m|-r-l(a:) (3.23) 

(|m| — r — 1)1 ' ' 



(cf. (EH) and 



Va,,. ■ ■ ■ Pa.,x (e-"'F|^|„,(x)) = e-'^'Uraix) (3.24) 



4. Non-null case 

4.1. Multiple Laguerre kernel. For the determinantal point processes (|1.4I) . the method of 
biorthogonal polynomials iF gives that the corresponding /c-point correlation function has the 
determinant form 

P{k){xu---,Xk) = det[/C(xj,j;j)]ij=i,...,fe (4.1) 

for a certain IC{x,y) referred to as the correlation kernel. The correlation kernel can always be 
written as a double sum over ii{x)r]j{y) weighted by a factor proportional to the inverse of the 
matrix of inner products [(Cij^j)]- We know from |1] that for the Laguerre case of ()1.4|1 this 
double sum can be written as a double contour integral 



ci-^> C{„i} 



1=1 ^ 



'b 



In the latter expression, a stands for a non-negative integer while C^ and Cr^i stand for 
non-intersecting counterclockwise contours that encircle b and —1, respectively. 
On the other hand, a determinantal point process 

1 ^ 
-7^Ww{xi)det\pi{xj)]ij=i^...^N<^ei[qi{xj)]i^j=i^...^N (4.3) 

in which {pj(x)}, {qj{x)} have the biorthogonality property, has for its correlation kernel the 
single sum 

N 

{w{x)w{y)f''^^Pi{x)qi{y) (4.4) 

j=i 

(cf. the first line in (|2.3|) ). Furthermore, the general theory of multiple orthogonal polynomials 
[TT] tells us that with Pi,qj as in 1)3. 7|) . (|3.8() can be summed according to a generalization of 
the Christoffel-Darboux formula in (j2.3|) . In the Laguerre case this summation has been made 
explicit in ^ . We revise this latter formula, and make explicit the form of (|4.4|) in the following. 
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i-1 

Proposition 4. Let Oj = (0, . . . , 0, 1, . . .). The correlation kernel \4-Sf^ has the summed series 
form 

(|m| + a — Ij! X — y \ ^—^ |m| ' / 

anrf i/ie single sum form 

Proof. First, we effectively multiply Eq. ()4.2p by (x — y)/{x — y) by writing 



/c^(x,y) 



27ri 7 27riu;-z Vl + w^y -y- V -^ - Oj / V dz dw ) ^ '' 

where use has been made of the relation between {6,} and {oj} noted in §3.1. Integrating by 
parts gives 



^Ux,y) 



2; X °/2 e(y-^)/2 



y J x-y 



dz f dw _,^^_^^y,^ f l + z ^^^'^^ / w-aA f N + a ^ 



-e 



27ri 7 27ri yl + tt;/ \^\z — aij \{l + z){l + w) ^ {z — ai){w — ai) 

which which from (|3.9|) is equivalent to the first equation of the proposition. 
Second, we apply Cauchy's integral theorem to Eq. (|4.2|) and obtain 



a/2 D 



UJ it ^ 2.i nf=i 



dz e^^^(l + z)l'"l+° 



-ic<-i} 



i 

{bi} 

D 



/dm py^ -X- 1 

27ri(l + ti;)IH+"ii^^~ ^'^ lo - z 

{-1} 



4-1} 



We can choose the contours C^. and C|_i} in a manner that guarantees \z — hi\/\w — 6j| < 1 
and |1 + w\/\l + z\ < 1. In such a situation, we have 

1 _y. {z-hiTH + wT 



w-z ^^ (w-bi)'^+^(l + z)'^+^' 
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and consequently, 



} 






C{-1} 



The integral in z around bi has no pole when j > rrii and so the sum over j can be truncated. 
We thus can write 

X;- (X y) = ixvr'' e-^^-^yy^vV I ^ e-(l + .)IH+-^--^ 






C{-i} 
Simple manipulations and recalling Hc{.9() complete the proof. D 

We now turn our attention to the correlation functions. It is known from jlj that they are 
given by (|4.1() with correlation kernel ()4.2() . To give further insight into this result, we will 
re-establish this fact using a formalism analogous to that revised in §2 for the null case. In 
particular we want to prove Eq. 1)2. 6|) remains valid in the non-null case. According to the 
discussion of §1 the perturbed Laguerre p.d.f. is given by 

1 ^ - 

P^(xi,...,x^) = ^n-"-"" n !^det[e--^];y^.^,, (4.5) 

JV i=l l<i<j<N ^ * 

where it is supposed that no eigenvalues in a coincide. When Oj = Oj for some i ^ j, the 
appropriate p.d.f. is obtained by applying L'Hospital's rule on P^{x). Further, we read off from 
(p^ that 

AT AT 

Z^ = (-l)^(^-i)/2iV! llia + i- 1)! n(l + a.)-^-"- (4-6) 

1=1 i=l 

Lemma 5. We have det [/C^(xj,Xj)][™i^ = |m|!i^^|(xi, . . . ,x^^^). 

Proof. Let us suppose ai < . . . < a^r where A^ = |m|; i.e., m = (1,...,1). Then Cauchy's 
theorem allows us to expand Eq. (|4.2|) . 

I — 1/7 7 — -L 

y) h nf=i («.-«.) J 2vriu;^+-ii/ 

feT^i WO} ■^■_^- 
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But we have 



dw e^^ y 



n-l 



2m w'^ (n-l)!' 



'{0} 



N 



N 



X{{w-l 









where a' = 1 + a = (1 + ai, . . . , 1 + cat) and where Cn (x) stands for the elementary symmetric 
function of degree n that does not contain the variable Xi dHl §1-3]. Combining the last equations 
we get 

N 

where 

^^'^■ = (a + j-l)! Vi(^)li(««-«fc) • 

We finally use 



(4.7) 



(4.8) 



and the general formula 



n(«;- - <) = det [(a',)^--i] .^. = det [e^;U^') 

i<j 



«J 



det 



'^Ck,ifk{xi)gi{xj) 



k,l 



det [/i(xj)]i_j- det [aj].^. det ^^(ajj)],^ 



«j 



to obtain the sought formula. 



D 



oo oo 



Proof. From Eqs (|4.7j) and (|4.7|) . we can write 



N 



N 



{a + j -1)1 



J, 7 = 1 



«J 



*'J=1 

which can be shown to be equal to N . Similarly 

oo oo 

dy/C^(x,y)/C^(y,z) = (xz)"/2e-(-+^)/2 

i,j,k,l Q 

\a/2 -{x+z)/2 ST^ „-aiXr^_ r^. _i-l("+J~l)' 



CXJ 

Y,e-'''''Ci,,Ck,iz'~^ j dye 



y-a^Vya+j-'^ 



{xzY'^e 



i,j,k,l 



(1 + afc)°+^' 



Direct manipulations then give 



7 ,CijCk,i-7^ 



a + j-l)\ 



(1 + afc)"+J' 



L^i.i 
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and the proof is completed. D 

Proposition 7. Let p"(xi, . . . , a;„) be the n point correlation function in the multiple Laguerre 
ensemble; that is, 

oo oo 

|m|l f f 

p'^{xi,...,Xn) := .1 I _ ^, / dXn+1--- / d2;|j„|P|^|(xi,...,X|m|). 



Then 

p^{xi, ...,xn) = det [ICZ,{xi,Xj)]l.^^ . 

Proof. On the one hand, we know from Lemma El that 

cx) oo 

p^(xi, ...,Xn) := _ / dxn+1 ■■■ j dx\^\ det [/C^(2;i,a:j)]J™i^ . 



On the other hand, Lemma IHl and the Laplace expansion of the determinant imply that 

oo 

j dxk det [/C^(xi,x,)]^_^.^, = (|m| + 1 - A:) det [/C^(xi,x,)];=;^, . 



Thus, starting with k = |m| and applying |m|— n— 1 times the last equation on det [/Cm(^i) ^j)]ij=i 
give the formula we wanted to prove. D 

4.2. Multiple Hermite kernel. The appropriate generalized Hermite kernel has been given 
in double contour integral form by Zinn- Justin I3H1. With our notation and normalization, it 
reads 

Cb -ioo *-^ 

where it is understood that z never crosses the lo's path. 

The following proposition will considerably simplify our analysis of the multiple Hermite 
kernel. It is proved by following the method exposed in Proposition [5J 



Proposition 8. Let X be a real number such that \\/2aX\ < a. Then, as a ^ oo, 

V2^ICZ,(a + V2^X,a + V2^Y) = ICn,(X,Y) + O ( ^] . 

Corollary 9. We have 



.1 g-(x2+y2)/2 / D 



^'"^^'^^ " 2H^(|m|-l)! x-y ( ^-(^)^-(y) " E ^^-+-»(^)^— >(y) 



and 

Proof. We first make the following change of variable in Proposition [^ x i— > a + \/2ax, a + 
v2ay, and hi i— > 6j/v2a. Then we use the asympotics of the multiple Laguerre functions given 
in Proposition [2 to write the scaled perturbed Laguerre kernel in terms of multiple Hermite 
functions. The proof ends with the comparison of the latter expression and the scaled kernel of 
Proposition IHl □ 
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In a similar way we can prove that the correlation functions in the perturbed Hermite ensemble 
are determinants of the multiple Hermite kernel. We recall that this ensemble is characterized 
by the p.d.f. 



N 



PNixu...,x^) = ^fle-^^ n ^^^det[e-»-^]f_. (4.10) 

The normalization constant Z^r is read off from (jl.Tj) . 






Corollary 10. Let 

oo oo 

|m|! f f 

Pn{xi,...,Xn) := .1 , _^s, / dXn+1--- / dx\^\P\^\{xi, . . . , X\^\). 

— oo — oo 

Pn{xi,. ..,Xn) = det [/Cm(a;i,a;i)]"j=i • 

Proof. The expressions of the Laguerre and Hermite p.d.f. (respectively given in Eqs 1)4. 5(1 and 
1)4.10(1 ) together with the Stirling approximation allow us to write 

lim^a2a)^/^P§{xi, . . . ,xn)[^^^^+,/2^^^ I = P^ixi, . . . ,xn). 

I bi^bi/V2a ) 

Hence the corollary is an immediate consequence of Propositions [T] and |S1 □ 

5. Quasi-null case 

In this section, we suppose that all but a fixed number of eigenvalues of the matrix A are 
zero. Precisely, 

a = (ai,... ,ar,0, ... ,0), 

b ={bi,...,bd,0), 

ui = {mi, . . . ,mfi,N — r). (5.1) 

Here d = D — 1 gives the number of distinct parameters in a that are not equal to zero. To 
make contact with the results of ^ and p2] reviewed in §1, we seek limiting expressions for the 
kernels, and thus for the correlations, when the size N of the random matrices becomes infinite. 

5.1. Perturbed Laguerre ensemble. The kernel is now given by 



w^ix) f dz f dw e-^^+y-' /l + zY+"fw^N-r^ yJ-a, 



where A^ = |m| and Wa{x) = x"e~^, as usual. It is convenient to introduce 



}C^ix,y) := .h^^}C^{x,y). (5.3) 

V Wa{x) 
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This gauge transformation does not affect the correlation functions because 



det 






-9[Xi,yjj 



N 



«J=1 



TV 



det [g{xi,yj)].. 



for all g{xi,yj) and f{x) / 0. 

A striking feature of the quasi-null soft edge scaled correlation (|1.12j) is that the correlation 
kernel consists of the null case Airy kernel (|1.14|) plus r perturbation terms. As our key result 
we now show that at the level of finite N, the multiple Laguerre kernel can be decomposed 
as a (null case) Laguerre kernel plus r perturbation terms. The latter are written in terms of 
incomplete multiple Laguerre functions of type I and II, 



and 



A(^)(x) 



A(^)(a 



-{0,ai,...,aj 



2vri^A.-nUi(^-«;^) 



dw e'^'^w^-'^ IY^-J^{w - ak) 



27ri 



(1 + t/;) 



N+a 



'{-1} 



respectively. Comparison with 1)3. 9|) shows that 

A»(x) = £;;r+i-^(x|b') 

when 

i 

W = {ai, . . . ,ai,0} and m' = {1, . . . , 1, A^ — r}. 
Similarly, we see from H3.1U() that 

(A + i-r)! 



A(^+i)(a 



(A + a-1)! 



X 



a+r—i—l—x r<a+r—i—l 

^ '"ml 



x\h') 



(5.4) 
(5.5) 

(5.6) 
(5.7) 
(5.8) 



As C^lx) and C'^{x) can be interpreted as anti-derivatives and derivations of Laguerre polyno- 
mials according to (|3.18|) and ()3.19|) . it follows from ()5.6() and (|5.8() that A(*^(x) and hS^{x) are 
also simply related to Laguerre polynomials. Explicitly, 



V ■ ■ -V 



a«(x)1=l^+_;+_\(x) 



and 



where in deriving the latter use has also been made of (|3.12jl . 

Proposition 11. Let K'^{x,y) be the gauged transformed Laguerre kernel; that is, 



f dz r dm p-^^+V'"' / 1 -\- z 

K^{x,y) = l^lC'^{x,y)= / ^ / |^^ (±±^ 

a^O / Z-Ki J ivri w — z V 1 + ^ 



for A^ = |m| and C\q^ n C|_i} 



Then 



C. 



N+a 



Z . 



{-1} 



^S.(x,y) = nt;(x,y) + ^A«(x)A«(y). 



(5.9) 



i=l 
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Proof. We first consider the relation 

1 -p-r 'u;-ai_ 1 , Y^ 1 lfc=i ■^ - «/= 

i=l 

For r = 1, this reads 



n^=^+E^f-^^^- ("») 



w-z^-z-ai w-z f^^Ul=iZ-ak 



1 w — ai 1 1 

+ 



w — z z — ai w — z z — ai 
which is trivially true. For r > 1, Eq. (|5.1Uj) is proved by induction; explicitly, 

r+l -, / r \ 

-■ ~ ^ — en \ z — a^+i + w — z 



_! TT W - aj _ 1 / T-r 

— y -^ ^ y — ft: ')n — r I -^ -L 



w — z-^-^z — Oi w — zX-^-^z — ai Z — Or+l 

1=1 \i=l / 

- ^ TT ^ ~"'i , [ TT w -aA I 

W — Z \^ Z — ai \ _ Z — Oi I Z — Or+l 



1 "■ 



Ul=i w-ak ^ Uk=i w-ak 



w — z 



r+l 



1 nUi z-ak Wk=i z-Ok 



as expected. 

We then substitute Eq. (|5.1U|) in the quasi- null kernel, given in Eq. (|5.2j) (see also Eq. (|5.3jl l 
This gives 

^„ , , , dz f dwe-'^'+y'" fl + z\^^'' ^w\N-r 



27ri ./ 27ri w — z \1 + w J \ z ■ 

dz e-^'^(l + z)^+° f dwey'^w^-''lli~JiW-ak 



^{-1} C{_i} 



+ E 



j=i 



27riz^-r]^*^^^(^_afc) J 27ri {1 + w)^+" 



Comparing the latter equation with Eqs 1)5. 41) . (|5.5|) finally completes the proof. D 

5.2. Perturbed Hermite kernel. When the set a of eigenvalues is quasi-null, the multiple 
Hermite kernel becomes 

ioO 9 / 9 ; 

J 27ri J Iiri w — z \ z J ^^ \z — a-i ) 

Note that the w's contour integral does note pass through the origin. As in the perturbed 
Laguerre ensemble, we introduce incomplete multiple Hermite functions. 



r«(^):=e--v2 / — — ^::ri^ — (5.12) 

''{0,ai,...,aj 



and 

ioo 



i-1 

r«(x) := e^'/2 /" ^e--W/4^A^-r TT^^ _ ^^)_ (5^^3) 

J 27ri -'--'- 



^==1 
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Comparison with (|l'{.2()|) and (|H.21j) shows 



f«fx) 



-x'^12 



(_2)A^+i-'-i(iV + i-r-l)! 



n^,{x\cl), 



( _-i \N+i~r 



IT 



(5.14) 



where b' and m' are the restricted sets of parameters given in (|5.7j) . Consequently, we read off 
from ()3.23|) . H3.24() that these special multiple Hermite polynomials are related to the classical 
Hermite polynomials by 



and 



f_l\N+r+i+l 



r«(a 






/-IT 



-.y2jj 



N-r[Xj 



The algebraic property given in Eq. (|5.1fl|) can be used to prove the following proposition. 
Proposition 12. Let Kj\f(x,y) be the usual Hermite kernel; that is, 



i^^(x,y) = lim^/C^(x,y) = e^'/2-V2 f ^ 

a— >0 J ZTTl ./ ZTTl 



Z J 



w — z 



Then 



/C„.(x,y)=K^_,(x,y) + j;r«(x)r«(y). 



1=1 



Before considering the large A^ limit of the kernels, let us link the perturbed Laguerre and 
Hermite ensembles by relating the incomplete multiple Laguerre and Hermite polynomials. 



Proposition 13. Set X a real number satisfying \\'2aX\ < a. Then, as a ^ cxd, 

(2a)(i+^-*-^)/2X»(„ + ^2^X) 

and 

(2a)(^+^-^)/2A«(a + ^/2^y) 



ai\—rai/y/2a 



^X-/2f{r)^X) + 0(-^ 



g-yV2r»(y) + 0('-^y 



a,;t— >ai/\/2a 

Proof. We essentially follow the steps given in the proof of Proposition |21 We skip the detail. D 
The latter result obviously complies with Proposition |B1 and the asymptotics 



l wa{a + y/2^X) ^ ^(y2_x2)/2 ^ q ( J_\ 

6. Kernels at the edges of the spectrum 

6.1. Hard edge. We know from (|1.14|) that the null-case hard edge kernel can be written in 
terms of Bessel functions Ja{x)- The latter satisfy 

dz e^(^-^"')/2 



Ja{x) 



27ri z' 



a+l 



J-aix) = (-l)"Ja(x), a £ 



(6.1) 



'{0} 
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The form (|5.9|) of the perturbed Laguerre kernel suggests that its hard edge scahng can be 
expressed in terms of the Bessel like functions (|1.20|) , (|1.2H) . Use of (|3.15j) and (|6.1|) shows these 
functions are simply related to the Bessel function by 

and 

The functions j'*^ and j'*^ are in fact incomplete versions of the multiple Bessel functions of 
type I and II, which we respectively define as 

_ _ /> J \m\+a-l ~x^z+l/{4z) 

JZi^) = JZi^H := {2xr I T^^TTT^ ZZ— (6-2) 

and 

jzix) = j^{x\u) := (2x)-" J ^- — ^i^;r.+r ^" • (6-3) 

C{o} 

Recall that D denotes, as in Eqs. (|3.1|) and (|3.2|) . the number of distinct parameters in h = i/"^. 
The link relation between the incomplete and the complete functions is obvious, 

J"it;+'"X^/^) = (2V^)°+'^+'"V«(x) and J°tr;(^/^) = (2V^)-"-'-+V«(x), (6.4) 
where {l'^} denotes a composition of length k in which all entries being equal to unity. 

Lemma 14. Suppose that Eq. (|5.1j) holds. Let hi = ANvi for i = 1, . . . ,d and /x = (m-i, . . . , rud)- 
Then, as N ^ oo, 



' J 2vri 


w^=i{z-v,r^ 


r (itye^'™- 


-l/(4»)j|^^(^_^.)n^. 


J 27ri 

C{o} 


y;|m|+a+l 



(4iV)-"£^(X/4iV) = (2VX)-"J-°(VX) + O (1/iV) , 
(4iV)-"£^(X/4iV) = (2^/X)-"J-°(Vx) + O (1/iV) . 
Similarly, for Oj = 4A^/ij (i = 1, . . . .,r), we have 

(4Ar)i-"-'-iA»(x/4iV) = jW(X) + O (1/iV) , 
(4iV)"+'^~*A»(X/4iV) = jW(X) + O (1/iV) . 

Proof. From the definition of >C^, given in Eq (|3.9j) . we have 

dz e-^^/4^fl + z)^+"-i 



CLiX/AN) 

'-{0,bi,...,bd} 

Then, changing z for 4A^tt; and 6j for ANvi, we get 

,nife 



£^(X/4iV) = (4iV)^ / o.. 



C{o,hi,...,h^} 



2^i ntii^-'^fc: 



But for iV ^ cx), 

(1 + l/4iV'u;)^+" = (1 + l/4A^u')^ + O (1/iV) = e^l^"" + O (1/iV) . 



22 PATRICK DESROSIERS AND PETER J. FORRESTER 

This immediately implies that 

~ /■ rim p-Xw+1/iw a+r-1 

''{0,hl,...,hrf} 

Comparing with Eq. H6.3() gives the desired result. The proof for the other functions is similar. 

D 

Theorem 15. Suppose that Eq. (|5.1() holds. Set Oj = ANhi for i = I, . . . ,q and Oj = {AN)'^hi 
for i = q + 1, . . . ,r. Following the notation of §5, let h := {hi, . . . , hq} = {y^\ • • • , fp''} =: f^ 
with q = |/x|. Then, as N ^ oo, 

4]v^S,(§,^)=C'^(^,n + o(l). 

where 

K};:^\X,Y) = (I) "'"' Kl^;^{X,Y) + J2J^\X)J^\Y) (6.5) 

n 7377 (6.6) 



27ri J 27ri w — z Kw/ ^_^ \ z — v,, 

^ ' '"{0} 



1 fY 



y" (jiX,Y) -2j2^^^J^'+.S^)J^-e.iVYyj . (6.7) 



2{X-Y)yX, ^ ^^^ 

with a' = a + r — q and 

J{X,Y) = J^'iVx) ((a' + q)J^'iVY) - VyJ^'+\Vy)) 

+ J^'iVY) ((a' + q)J^{^) - VX^;'-i(VX)) . (6.8) 

Proof. Let us begin with Eq. (|6.5() . We know from the previous lemma that 

p)(X/4iV)AW(y/4iV) = (4iV)jW(X) j(*)(y) + 0(1) when i=l,...,q. 
For i > q, Lemma [TBI and Eqs (fTTT^ - lfT^ imply 

A{i+9)(X/47V)A(^+5)(y/4iV) = jW(X)j('^)(y) + Or^") when i = l,...,r-q. 

Furthermore, we have from Eq. I\b.'d\i and §2 (see also |151 §4]) 

1 / Y V \ / Y\ ("+'')/2 / 1 



where the hard edge kernel is defined in Eq. H1.19() . We finally obtain Eq. 1)6. 5p by using 
Proposition 111! 

We now turn our attention to Eq. 1)6. 6p . We first claim that the hard edge kernel has the 
double integral representation 

K'ri.,y) = (i) I^J '£'- ^^ (i) . 

Cloy C^""^ 
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This is proved by effectively multiplying the previous equation by {x — y) / {x — y) and integrating 
by parts; that is, 



X ^ y \y / J ^""i J ^^1 w — z \wy \ oz aw J 



g xz+yw\ 



i/2 



1 (x\ ' f dz f c?u)^i/(4^)_i/(4^)^_^^+j^„ /z\" / a z + w 



^ ^ y \y J J 27ri J 27ri Vw/ V^zw 42^1/;^ 



'{0} 



^{^} 



'{0} 

1 



2aJa{^/x)Jaiy/y) - y/xJa-l{Vx)Jaiy/y) - JaiVx)y/yJa+liVy)) 



2{x-v) 

which, by virtue of ■v/xJ4(\/x) = ^/xJa-l{^/x) — aJa{^/x) and y/xJa+i{y/x) = —^/xJa-l{^/x) + 
2aJa{\/x), turns out to be equivalent to Eq. p.l9j) . Then we use the integral representations 
of the incomplete multiple Bessel functions, which are given in Eqs (|1.2()j) ~ (|1.2H) . to rewrite 
Eq. (iniSl) as 



" ' J 27ri J 2m 



^{0}uh c\l\ 



\w) [ w-z + £- JJ- [z-hj 



where the Cauchy theorem has been applied for deforming the contours C/q} and CsQ^h-^^j^.y into 
C{o}uh- Use of Eq. (|5.1U|) then finishes the proof of formula Eq. (|6.6|) . 

We show Eq. (|6.7|) with the same method than the one explained in the proof of Proposition 
H we multiply Eq. (jHSl) by (X - y)-i(X - Y) = -{X - Y)-^e^''-^'"{d^ + du,)e-^''+^'"; we 
integrate by parts; and we compare the result with the definition of the multiple Bessel functions, 
Eqs (jH^-dlll. □ 

6.2. Soft edge. Working directly with the double contour form ()5.'2() of the correlation kernel 
Baik, Ben Arous and Peche j^ have shown that the soft edge scaling (|l.in|) of perturbed Laguerre 
kernel can be written in terms of multiple Airy functions according to the results H1.1U() - ()1.16() . 
From the contour integral representation of the Airy function, 

Ai(x) := / ^e-"^+^'/3, (6.9) 

J 2-711 

A 
and p.lSf) we see that the multiple functions ()1.14|) . (jl.lSp are related to the Airy function by 



Ai{x)=V.s,---V^s, 
and 



Ai^'Vxl 



According to §6.1, the functions Ai'*-* and Ai can be completed in order to get the multiple 



M('Hx) = {-iyVs,---V,^_AAi{x)]. 

ctio 
Airy functions of types I and II, 

Mm{x) = MmixH := f ^e--^+-'/'f[{v - a,)— ^ (6.10) 

^ 1 — 1 



D 



A. ^=1 
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Figure 2. Contours A and A 



{si,...,Si}- 




and 



Ai^ix) = Alexia) := (-1)1-1 / -^^e---+-'/-^l[(v + ak 

One can easily relate the two sets of multiple Airy functions, follows: 

Ais^ii^{x) = Ai (x) and Ais^ii-i^{x) = — Ai*-*^(x). 



lOlfc 



(6.11) 



(6.12) 



Lemma 16. Suppose that Eq. (f^ holds. Let A = AN and B = 2(2iV)^/3. Let also bi 
— 1/2 + ai/B for i = 1, . . . ,d, // = (mi, . . . ,md) and \fj,\ = r. Then, as N ^ oo, 



B'-'C'^iA + BX) 



(-1) 



Af+r+l 



2°+' 



--^Ai^^iX) + O 



Ari/3 



(-1) 



N+r 



1 



B^^^Ci^iA + BX) = ^-^^^.^^Ai^iX) + O ^^^/3^ . 
Similarly, when Oj = —1/2 + Si/B for i = 1, . . . ,r, we have 

^i-^g-(A+BX)/2;^«(^ + BX) = {-lf+'-2-''-^Ai'\x) + O 



1 



Qi^iA+BY)/2^ii)^^ + BY) = (-l)^+^^2"+^^Ai«(y) + O (^) • 
Proof. The proof is based on the saddle point method. For instance, consider 

A('){A + BX):= I 

C{0,ai,...,ar} 



|^e^/W5(.) 



with 



-BXz^r 



/(z) = -4z + ln(l +z) -Inz, g{z) 



z'-{l + zY 



Y{\=i{z-ak) 
The function / has a simple saddle point of order two at zq = —1/2; that is, 

/'(^o) = = /"(zo) r(zo) / 0. 
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Since /'"(zq) = 32, we have that {z — zo)f"'{zQ) is minimum when arg(z — zq) = — 7r/3, vr/S, tt. 
We choose ±7r/3. Doing the change of variable v = B{z — zq), we understand that the steepest 
descent contour in the complex v-plane starts at ooe'^™''^ crosses the origin and ends at ooe''^'^. 
Moreover, 

Nf{z) = 2N - mN + v^ /?, + O (n^^/^^ , 
while 

i 

g{z) = (_i)'^2-"-"5^e-^^+^^/2 ^{v- Sk)-^ (l + O (iV-i/3 

fe=l 

if Sfc = B{ak + 1/2). Therefore, 



27ri 

-4{si,...,sJ 



fc=l 



This proves the asymptotics of A'*^ . The asymptotic behavior of the other functions is proved 
similarly. D 

Theorem 17. Suppose that Eq. (|5.H) holds. Set ai = —1/2+Si/ B fori = 1, . . . ,q and a-i » —1/2 
for i = q + 1, . . . ,r. Following the notation of §5, let s := {si, . . . , Sq] = {o"^% . . . , CTp'} =: o■^^ 
with q = \fj,\. Moreover, let A' be a contour going from ooe'^^'^'^ to ooe~^"^'^ and not intersecting 
Act- Then, as N ^ cx), 

B)C^ {A + BX,A + BY) = lCf\X, Y) + (1) , 



where 



-Xi) 



1C°^\X,Y) = i^^°f*(X,y) + ^ Ai^ \x)Ai»(y) (6.13) 



i=l 



du I dve e -Q^!_^ ^q.U) 



27ri J 27ri v — u ^J- \u — a, 

Aa- A' *-^ 



^^\y) (^(^'^) + ^/X.^V+e,(X)^V_e,(l^) ) (6.15) 



i=l 



with 



AiX,Y)=Ai,iX)'-^!^-Ai,iY)'^^'^^ 



dY '^' ' dX 

Proof. We first consider Eq. (|6.13jl . The previous lemma implies 



k^^{A^BX)k^\A^BY) = B-^Ki^\x)Ki^^i^)^o{A^ when i = l,...,g. 
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For i > q, the supposition Oj » —1/2 is equivalent to Oj = —1/2 + Si/B together with Sj = Bsi 
for some positive real number Sj. Then, Lemma [THl and Eqs (|1.14|) - (|1.15j) imply 

X(i+i) (A + BX) A(^+'?) (yl + BY) 

= i?-2Ai^'\x)Ai(«)(y) + o(^) =o(^) when i = l,...,r-g. 
Moreover, it is well known (see §2) that 

BK%{A + BX,A + BY) = K'°^\X, Y) + O (j^) ■ 

From Proposition ^2 we conclude that 

1 



BIC^{A + BX, A + BY) = BIC^{A + BX, A + BY) + O 



Ari/3 



BK^t'^iA + BX,A + BY) + bJ2 A^*) (a;) A« (y) + O (j;^) 

i=l 

K^°'\X,Y) + ^ Ai^^X)Ai»(y) + O (j^) 
1=1 



as desired. 

We now prove Eq. (|6.14|) . By using integration by parts, it is a simple exercise to show that 
the soft edge kernel defined in Eq. (|1.14|) can be written as a double integral, 

r rill r Hi) p-Xu+u^ /■i Yv-v'-^ /d, 
K^°'\X,Y)= /^ / ^^^ 



27ri J 27ri v — u 

A A' 



Also, the integral representation of the incomplete multiple Airy functions, given in Eqs (|1.15|) - 
HHni), leads to 

^ f ^^ ^-XuW/3^Yv-v^/3 Ul~=liv - Sk) 



S'",x)Ai<..,K). I ^J 



A{si,...,si} A 



2vri7 2vri UUii^ - Sk) 



Cauchy's theorem allows us to deform the contour A{si,...,Si} iiito A{si,...,sr} — -^{o-i,...,CTd}- 
Eq. (|6.13|) then gives 



9 T-ri-l 



lC'::'\x,Y)= I — /■i!ie-^«+«V3gy.-.V3(^ + y- 

Eq. ()6.14|) is then obtained by exploiting the decomposition (J5.1U|) . 

We finally get Eq. (|6.15;) by multiplying Eq. (|6.14|1 by {X - Y)/{X - Y), which gives 



X — Y J 27ri J 27ri v — u ^^ \u — aij \ du dv 

A„ A' '=^ 

by integrating by parts, and by comparing the result with the definitions of the multiple Airy 
functions given in Eqs (|6.1fl|) - (|6.11j) . D 
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The ease at which the hniiting kernel (|6.1H)) can be deduced from Eq. (|5.9j) contrasts with the 
difficulty encountered in j^ in computing the same limiting kernel from a double contour integral 
representation equivalent to (|5.2j) . The reason is that in the double contour the saddle point 
is such that w = z in the integrand, but then the denominator vanishes. By having reduced 
the double contour to a sum of products of single integrals as in Eq. (|5.9j) , the complication is 
avoided. 

We finish this subsection by showing that the kernel of the perturbed Hermite ensemble can 
also be mapped to the soft edge kernel /C^^*. By doing this, we generalize a result previously 
obtained by Peche EH] (i.e., Eq. KT^ for si = . . . = s^ = 0). 

Proposition 18. Let K,m{x,y) be the kernel given in Proposition^^ Let also A = \/2N and 
B = l/\/2iVV3. Suppose ak = \f2N{-l + Sk/N^'^) and Sk > for k = l,...,g. Suppose 
moreover Ok ^ —A for k = q + 1, . . . ,r. Then, as N ^ oo, 

q 

BIC^{A + BX,A + BY) = K"°f*(X, Y) + Y, Ai (X) Ai^^) {¥) + O (n'^/^ 

i=l 

Proof. This formula is a straightforward consequence of Proposition ^J and the following as- 
ymptotic relations, 

f «(A + BX) = i-lf^^^^^^^^^^^\x) (l + O (iV-V3)) , 

AN+i-r -N/2 

r»(A + BX) = (-1)^+'+^- -1 Ai«(X) (^1 + (^N-^/'^ 

These equations can be obtained from the integral representation of the incomplete Hermite 
functions by using the steepest descent method (see for instance Lemma 115(1 . We skip the 
detail. D 

6.3. From hard edge to soft edge. It is known (see for instance |151 §4]) that we can go from 
the hard edge to the soft by rescalling the spectral variables and taking the limit a — > oo: 

2a(a/2)i/3^^ard(^^ y^ ^ i^^°ft(^, j]) + O (a~^^^) (6.16) 



\i X = a^ — {20^)"^'^^ and Y = a^ — (2a^)^' ^77. This formula is a consequence of the asymptotic 
relation between the Bessel and the Airy functions I28j . 



(I) UVx)=Km + 0{a 



In the subsection, we show that the mapping between the hard and the soft edges is preserved 
in the perturbed Laguerre ensemble. 

Lemma 19. Let X = o? - {2c?Yl'^i and avi = 1/2 - ai/(4a)^/3 with ctj > for i = 1,. . . ,d. 
Then, as a ^ 00, 

Q,N (1-|m|)/3 ~ . — I I ~ /I 

I) J^{VX\,.) = {-1)M Ai^i^la) + 0[Jj^ 

q,n(1+|m|)/3 ^- 11 /I 

j;:{^\u) = {-l)MAi^{^\a) + 0' 



2) -M^v^M"; V -J ^^"Mvsi"; ^Val/3 
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Proof. The proof relies on saddle point method for both formulas. Here we only prove the first 
asymptotic development. We have from Eq. (|6.2|) 



j;(^/X|r.) = (2a)°(l-(a/2)V3^/ar J 



C{o}ui> 



^^ ^\fj.\+a-l^-a^z+{2a'^)^/'^liz+l/{4z) 
2^ ULA^^^' 



By making the change z i-^ z/a and by considering (1 — (a/2)-'^' ^^/a)" = e '"'^' ^ + (a ^'^^ , 
we get 



^{0,ai^l,...,aiyj^} 

where 



g(4a)l/3^2^|^|_l 

f{z) = —z + l/Az + lnz and ^(z) 



The function f{z) has a second order saddle point at zq = 1/2. We now suppose that aui < 
1/2 for all i. Following the method used in proving Proposition I16L we deform the closed 
contour C^Q^aui,...,aiya} ^^ such a way that it reaches the steepest descent contour, which is a path 
approaching zq from 9z < with an angle of — 27r/3 and leaving zq for Qz > with an angle of 
27r/3. Setting 

z = -v/{Aay/^ + 1/2 and aui = 1/2 - ai/{4a)^/^ , 

we find that 



a 



f{z) = -2a In 2 + ^^3 + O (a^^/^ 



(_l)lMl2alMl/3 e-«.e(a/2)V35 
Returning to the integral representation of the type I multiple Bessel function, we get 

Comparison with Eq. (|6.10() shows that this is the expected result. D 

The above result can be applied to the incomplete Bessel and Airy functions. For instance, 
one can show from Eqs ()6.4() and (|6.12|) that 

2a(«/2)V3 (^^y ji^{X)J^\Y) = Ai^^e)Ai«(r?) + O (-1^) . 

When substituted into Theorem 1151 this formula together with Theorem El imply the following 
asymptotics. 

Proposition 20. Let X = o? - {2a^f/^i, Y = a^ - {2a^f/^r] and avi = 1/2 - ai/iAaf/'^ for 
ai > 0. Then, as a ^ cxd, 

2a(a/2)V3/ch-d(j^^y|^) ^^C^oft^^^^l^) ^O / -1/3 
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7. Concluding remarks 

— {«) 
7.1. The case r = 1 at the soft edge. Noting from (|1.15() that Ai (x) ^ as x — > — oo, we 



see from (|3.15() that 



X 



Ai<"m./e..<-.A.W... 



Substituting in (|1.12|) shows that for r = 1, si = d, K,^ {X,Y) can be written exphcitly in 
terms of Airy functions 



X 

IC'°^\X, Y) = K'°^\X, Y) + Ai(y) / e^(^-^) Ai(x) dx. (7.1) 

— oo 

In fact 1)7. 1|) has previously been derived as the soft edge scaled correlation kernel for random 
matrices closely related to p.8|) J16| I17|. These are (A^ + 1) x A^ complex Gaussian matrices X 
with the entries of the first A^ rows independent complex Gaussians with mean and variance 
1/2, while the entries of the final row are distributed according to this same complex Gaussian 
multiplied by ^/c. The quantities studied were the eigenvalues of X^X. With c = 1/b, such 
random matrices differ from those specified by (|1.8|1 . and the paragraph including (|1.1()|) with 
r = 1, only in that the final row of X has variance l/2yb instead of the final column. 

Results of 16, Prop. 14] give that with the largest n eigenvalues {xj}j=i „ of X^X scaled 
by 

Xj ^AN + 2{2Nfl'^Xj (7.2) 

and with 

l + A d 

the corresponding n-point correlation function is given by 

det[iC'''\Xj,Xk)\j,k=i,...,n (7.4) 

where 

Y 

^soft(^^y^^_gd(x-y)^j^g-dx f e'^^K'°^\X,v)dvy (7.5) 

~oo 

The scalings ()7.2|) , ()7.3|) are identical to (jl.lUj) . Since we have shown that the soft edge scaling 
for the perturbed Laguerre ensemble is independent of the parameter a, one might suspect that 
it is independent of variance being altered along a row instead of down a column. Indeed use of 
the integral form of K^°^^ (|1.14|) . 



oo 

K''°^\X, Y)= I Ai(A + t)Ai(y + t) dt, 



and use of integration by parts in (|7.5j) reproduces (|7.1j) . but with A ^^ y, an operation which 
leaves unaltered (|7.4|1 . 



oX V / a=0 / 
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7.2. An open problem. For the same class of random matrices which gave rise to (|7.5)) . the 
n-point correlation function for the smallest n eigenvalues {xj}j=i,...,n, scaled by 

Xj ^ Xj/4N, b = 2Nw 

was computed in ^1 Prop. 14] as being equal to 

det[/C '''' {Xj,Xk)]j^k=i,...,n 
where 

Y 

d_ 
dXV 



We have not been able to express 1)7. 6|) in the form of H1.18() with r = 1. Thus unlike the 
situation at the soft edge, the eigenvalues at the hard edge distinguish the variance of the final 
row of X being l/vb rather than the variance of the final column. 

This is not surprising when one recalls that for matrices W^^BW with W as in H1.8|) and B 
an n X n positive definite matrix, the eigenvalues p.d.f. is given by |in| I.SDj 

\\j=lJ- i-i-j<k\'^i ~ "k) L k=l,...,N-n k=l,...,N^ 

and is thus different in the neighbourhood of the smallest eigenvalues to (|1.9|) . A challenge for 
future studies is to compute the n-point correlation for the hard edge scaling of the smallest 
eigenvalues specified by (|7.7|) with q = n — A G Z"*" general. 
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